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THE SCATTERING OF POLARIZED LIGHT BY 
POLYDISPERSE SYSTEMS OF IRREGULAR PARTICLES 

By Alfred C. Holland 
Electronics Research Center 


SUMMARY 

The elements of the Mueller matrix for polydisperse systems 
of irregular, randomly oriented particles have been measured in 
absolute terms as a function of scattering angle for one wave- 
length. These results have been compared to the matrix elements 
that were calculated for assemblies of spherical particles that 
fit the same particle size distribution function and have the 
same (real) refractive index. Correlations between the measured 
and calculated matrix elements are discussed. 

I.- INTRODUCTION 

The atmospheric aerosols exert a profound influence on the 
optical properties of the earth's atmosphere. Despite their low 
concentration relative to the molecular constituents, the aerosol's 
effect on atmospheric scattering can never be completely ignored, 
particularly when polarization phenomena are considered. The 
atmospheric aerosol is composed of particles of many different 
shapes, sizes, and materials. Most of the aerosol particles 
(dusts, smokes, snow, and ice crystals) are known to have non- 
spherical shapes. 

Before a complete description of light scattering by the 
atmosphere can be made, one must know the detailed scattering 
behavior of such nonspherical particles: specifically, how the 

scattering matrix of the atmospheric aerosol depends on the 
particle composition, size distribution, shape, and, finally, on 
the wavelength of the incident light. However, except for 
spherical particles, and within certain limits, ellipsoids and 
cylinders , no rigorous theoretical description of particulate 
scattering is known. 

The practical importance of the problem can hardly be over- 
estimated. Beyond the case of atmospheric optics and the effect 
of the atmosphere on solar insolation and, hence, climatology, 
there is the possibility of using light beams to probe the 
atmosphere. Determining the properties of an unknown medium by 
examining the behavior of its scattering matrix with scattering 
angle and wavelength has many potential applications . 

At the present time, most investigators (for lack of any 
better practical model) assume that their scattering particles 



are spheres. The scattering properties of spheres can be 
calculated exactly with a large digital computer. 

Deirmendjian (ref. 1) for example, has calculated scattering 
matrix elements for models of haze and fog and found his results 
to be in excellent agreement with the measurements of Pritchard 
and Elliott on natural fog (ref. 2). 

There has been some experimental evidence that the scattering 
behavior of assemblies of irregular but randomly-oriented 
particles is similar to the scattering behavior of assemblies 
of spherical particles. Hodkinson (ref. 3) showed that the angular 
distribution of light scattered in the forward direction by 
samples of quartz, flint, and diamond dusts was similar to the 
light distribution predicted on the basis of geometric optics. 
Holland and Draper (ref. 4) showed reasonably good agreement 
between the unpolarized scattering coefficient for polydisperse 
samples of talc and that calculated for spherical particles 
fitting the same size distribution. 

The experiments reported here have two main objectives. 

The first objective is to investigate possible correlations 
between the scattering matrix elements of irregular (but randomly 
oriented) particles and the matrix elements calculated for 
spherical particles having the same size distribution and 
refractive index. The second objective is to compile sufficient 
data on the scattering properties of irregular particles of known 
size distribution and physical properties to aid in the develop- 
ment or testing of the necessary theoretical models. 

Considering the second objective, this work can be viewed 
as a preliminary step in a program aimed at developing methods 
for describing the scattering of light by real media. 

II.- THEORETICAL CONSIDERATIONS 

The scattering of a polarized beam of light by an optically 
thin volume element of a scattering medium can be described 
efficiently by the Mueller calculus. The scattering event is 
described as a linear transformation; the Stokes vector of the 
incident light beam is transformed into the Stokes vector of the 
scattered beam by the Mueller or scattering matrix: 

{h} = oonst [ s i j ]{ r j } • 

The Stokes vectors of both incident and scattered beams depend 
only on the characteristics of the two beams, while the scattering 
matrix depends only on the characteristics of the scattering 
medium. 
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The Stokes parameters and the Mueller matrix have been 
discussed fully by many authors including Perrin (ref. 5) , Van 
De Hulst (ref. 6) , Rozenberg (ref. 7) , and Parke (ref. 8) and 
will not be described in detail here. In brief, if one chooses 
as the reference plane the scattering plane, defined by the 
directions of the incident and the scattered beams, then the 
complex electric vector of either beam can be decomposed into 
components perpendicular and parallel to the reference plane. 
With this convention, the Stokes parameters are defined as: 

< ' * * \ 

E.E„ + E E > 
s Z Z r r / 

Q - < E * E <1 - E r E r> 

u = <KK E ri>= ( E l E r +E r E l) 

V = <-2Im jE^O = (i< E A +E rV) 

The subscripts £ and r denote components parallel and perpendic- 
ular to the reference plane; the brackets denote time averages 
of the noted quantities, and the asterisks denote complex con- 
jugates . 

The Stokes parameters all have the dimensions of intensity 
and satisfy the following inequality: 

o 2 2 2 

I* >_ Q + U + V . 

The equality holds only for completely polarized light. It can 
be shown (ref. 9) that this characterization of a light beam is 
unique and that two beams of light with identical Stokes para- 
meters cannot be distinguished by optical analysis. 

The four Stokes parameters can be considered as components 
of a four-vector and a scattering event described as a matrix 
multiplication : 
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In general, the scattering matrix has sixteen independent 
coefficients; however, Perrin has shown that if the scattering 



medium is isotropic and free of optical activity the scattering 
matrix has only six independent non-zero coefficients: 



The form of the matrix for different types of media has been 
discussed in detail by Van De Hulst (ref. 6 ). Rozenberg (ref. 7) 
points out that in general the matrix elements depend not only 
on the material properties of the medium but on the scattering 
angle 0 measured in the scattering plane and on <f> the azimuth of 
the scattering plane itself from some reference direction. Only 
in the case of an isotropic medium can one define matrix elements 
that are independent of the azimuth angle <j> and that depend only 
on the usual scattering angle. 


For a single homogeneous isotropic sphere we have the 
additional simplification Sjq = S 22 and S 33 = S 44 . There are, 
then, only four independent non-zero matrix elements. These 
four coefficients can be expressed in terms of the intensity 
functions ij(m,x,0), (j=l,2,3,4) and the scattering efficiency 

Q s (m,x) obtained from the Mie theory: 
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S 

S 
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1:L (m,x, 0 ) 

= S 22 (m,x , 0 ) 

12 (m,x, 0 ) 

= S 21 (m,x,0) 

33 (m,x, 0 ) 

= S 44 (m ,x , 0 ) 

43 = “ S 34 

= i 4 (m,x, 0 ) 


^■|i 2 (m,x, 6 ) 
^-|i 2 (m ,x , 0 ) 
i 3 (m,x, 0 ) 


+ i 1 (m,x, 6 ) 
- i^ (in ,x , 0 ) 


The matrix elements defined above hold for a single spherical 
particle of refractive index m and Mie size parameter x = ttD/A . 
The parameter D is the diameter of the particle and A is the 
wavelength of the incident light. 

For a polydisperse collection of spherical particles whose 
diameters are distributed according to a frequency function 
f (D) dD (particles/cm 3 -cm) , one can define comparable volume 
intensity functions as: 
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Pj (xn,x o , 6) 
4tt 



(m,x) 



f (x) ij (m,x , 6) dx 


( j = l ,2,3,4) 


The volume scattering coefficient is defined in terms of the Mie 
scattering efficiency: 

~ X M 

2. 


6 S (m,x) 


= — f 

^ J 


x^Q (m,x) f (x) dx 
s 


In the expression above we have used x & = TTa/X , x M = ttD m /A; 
k = 2ir/X is the optical wave number or propagation constant. The 
quantity x Q is the Mie size parameter corresponding to the median 
particle diameter of the size distribution; x Q = ttD 0 /X . With 
these definitions, the scattering matrix for a volume element of 
homogeneous, isotropic spheres that fit the distribution function 
is : 
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The matrix elements are normalized so that 
J 2] 4ti 

The particular particle samples were found to fit a general 
lognormal distribution function: 


x o' 6) 
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4tt 


do) = 1.0 


f (D)dD 


N6 

(D-a) 




particles 
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with the constraint a <_ D <_ °°. The pertinent parameters are 
the median particle diameter D 0 = 0.45 microns, a = 0.05 microns 
and 6 = 0.896. 
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For our purposes, it is convenient to express the distri- 
bution function in terms of mass concentration (gms/cm 3 ) rather 
than particle concentration (particles/cm 3 ) . For spherical 
particles, the particle concentration N is related to the 
particle mass concentration M through the equation 


tt P E (a , 6 ) * 

The parameter p is the mass density of the particle material 
and the function E(a,6) is: 

E (a , 6 ) = a 3 40a 2 (D -a) exp j-^l+3a (D -a) 2 exp j— ^-1+ (D -a) 3 exp/-A T 

° '46 ' '6 } '4S 2 


The volume intensity functions for our spherical S . 0 0 model are 
then: 1 2 
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and the volume scattering coefficient is: 
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The mass scattering coefficient is defined as: k s (m,x 0 ) = 

B s (m,x 0 )/M with dimensions of cm2/gm. 


The preceding equations were programmed for the IBM 7094 
computer in Fortran IV language. The single particle properties 
were computed for values of x = 0(0.2)120 for scattering angles 
0 = 0° (1°) 10° (2°) 30° (5°) 160° (2°) 180°. These single particle 
properties were then integrated over the size distribution function 
using a step size of Ax = 0.2. Convergence of the integrals was 
checked by the program. 


III.- SCATTERING MEASUREMENTS 


All of the quantities that are necessary to determine the 
scattering matrix elements of the particle samples can be measured 
with a recording Polar Nephelometer . The instrument used in 
these measurements is a model PN-2 recording Polar Nephelometer 
manufactured by Spectrametrics Inc. of Burlington, Massachusetts. 
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Basically, the instrument consists of a source telescope that is 
fixed relative to a circular mounting base; a scanning telescope 
that is mounted on an arm pivotted at the center of the base; and 
a fixed telescope mounted so that it monitors light scattered at 
an angle of 11 degrees. The scanning telescope covers scattering 
angles from 18 degrees to 166 degrees. Light traps are mounted 
opposite each telescope to provide a uniform black background. 
Overall diameter of the instrument is 8 feet. The essential 
features of the instrument are shown in Figs. 1 and 2. 

The detectors used are end-on photomultipliers with S-20 
photocathodes. Polarizing filters and spectral filters can be 
mounted on both of the receivers and on the source. The light 
source is chopped at 90 cps , and synchronous detectors and 
adjustable time constant integrators are used in both receivers. 
The optical systems of both the source and the receiver telescopes 
are identical. The source beam and the field of view of both 
receivers are rectangular, roughly 2 inches high and one-half 
inch wide. The intersection of the source beam and the receiver 
field of view defines the common volume. As the scanning arm 
moves, the volume changes size and a correction for this effect 
is made during data reduction. The time required to scan through 
the complete range of scattering angles can be varied from 1 to 
over 20 minutes. 

Calibration of the Polar Nephelometer System 

The system was calibrated in absolute quantities by comparing 
the intensity of the light scattered from the particles in the 
common volume with the intensity reflected or transmitted from a 
calibration standard. The method is identical to that used by 
Pritchard and Elliot (ref. 2). Since a complete discussion and 
proof of the method are given in this reference, only a qualita- 
tive description is given here. Briefly, a calibration standard 
made of thin plastic sheet of measured reflectance and trans- 
mittance is placed in the common volume perpendicular to the 
scattering plane. The calibration screen is moved slowly in a 
direction perpendicular to the screen surface so that each pair 
of rays common to the source and receiver beams is intercepted 
during the screen's travel. The intensity of the light reflected 
by the screen is recorded as it moves , and the integral of the 
signal over the distance travelled is recorded. This integral is 
compared to the intensity of light scattered by the particle 
samples. The reflection and transmission matrices of the 
calibration screen are measured and compared to an absolute 
standard; a glass plate covered with magnesium oxide smoke 
particles. The reflection matrix of the MgO standard has been 
measured by Pritchard and Elliot (ref. 2) . 




Figure 1.- Schematic of polar nephelometer 





SCATTERING 

ANGLE 


( 0 ) 


Figure 2.- Block diagram of signal processing system 
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The calibration was carried out for every combination of 
scattering angle, polarizing and analyzing filters, and color 
filters used in the measurements. The procedure yields a 
calibration constant K(P,A,0;A) such that the volume scattering 
coefficient for the sample volume is: 

B (P,A,0;A) = K (P , A , 0 ; A ) X W (P ,A, 0 ; A ) 

s SC 

where W sc (P,A,0;A) is the output signal from the scanning photo- 
meter channel for source polarizer, P; receiver analyzer A at 
scattering angle 0; and effective wavelength A. The effective 
wavelength is the center wavelength of the combination of source 
spectrum, filter transmittance and photocathode response. 

Particle Dispersion System 

The particles are carried through the common volume by a 
jet of dry nitrogen gas. The jet axis coincides with the axis 
of rotation of the scanning arm. An exhaust system scavenges 
the spent particles through the hollow hub of the scanning arm. 

The particles, because of their small sizes, take up the motion 
of the gas stream and the jet can be described as a two-phase 
gas jet. The mass fraction of particles in the stream is small 
enough so that their presence does not effect the flow properties 
materially . 

The particle density in the common volume was determined by 
surveying the velocity distribution in the jet at a sufficient 
number of points along the jet axis and radially from the axis. 

The velocity measurements were made both with and without particles 
in the stream. No essential velocity differences were noted. 

Then, associating the particle mass concentration in the jet with 
ordinary gas mixing rates and employing subsonic, two-phase, flow 
theory (ref. 10), the particle density distribution throughout 
the jet was calculated. A second method was used to substantiate 
the first calibration. Photographs of the jet illuminated by 
parallel light were taken and processed with a Joyce-Liebl 
microdensitometer. The density distribution in the jet was then 
determined using a standard Abel inversion (ref. 11). The 
relative density distributions obtained were identical to those 
given by the velocity measurements. With the jet calibration 
completed, the mean particle mass concentration in the common 
volume can be expressed in terms of the axial velocity of the 
jet at the center of the common volume and mass flow rate in the 
jet. Since the size of the common volume changes with scanning 
angle, the mean density of the particles in the volume changes 
and a correction must be applied during the data reduction. 
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Particle Characteristics 


For these experiments we used the Pennsylvania Glass Sand 
Company's "5 micron minusil" particles. According to the 
manufacturer's specifications, the material is a high-purity 
crystalline silica (99.9% SiC>2) with a specific gravity of 2.65 
and a refractive index of 1.55 for visible wavelengths. The 
particles have a flat plate-like shape with irregular edges. 

Some typical particles are shown in Fig. 3. 

The particle size distribution was determined by sampling 
the air-particle jet at the common volume. Electron micrographs 
were prepared from the samples and the size distribution was 
measured with a Zeiss Particle Size Analyzer, which permits 
measurement of particle size frequency and cumulative size 
distribution. Both functions were measured for our samples using 
well over 6,000 particles for each function. The data were 
found to fit a general lognormal function. The particle frequency 
function and cumulative distribution function are shown in Figs. 

4 and 5 . 


IV.- SIGNAL PROCESSING AND DATA REDUCTION 

The jet containing the particles is turbulent, and fluctu- 
ations in the instantaneous particle concentration result in 
fluctuations in the recorded signal. To eliminate these fluctu- 
ations, the fixed photometer was used to normalize the scan 
photometer signal. A block diagram of the signal processing 
system is shown in Fig. 2. The signal from the scanning photo- 
meter passes directly into a synchronous detecter or lock-in 
amplifier set at the chopping frequency. The signal is then 
amplified and scaled. A similar system processes the signal from 
the monitor photometer. Both signals are then fed into an analog 
divider whose output is a scaled quotient: 


const x W sc (P,A,0)/W m (P,A, 8 m ) . 

The output signal from the divider is displayed on one axis of 
an x-y recorder against scan angle. The scan angle signal is 
provided by a ten-turn, high-precision potentiometer driven by 
the scan arm. In addition, the monitor signal is integrated 
over the running time of the jet, and the resulting integral and 
jet running time are recorded. 

Since the monitor is fixed at a scattering angle of 11 
degrees, fluctuations in its signal are only due to jet density 
fluctuations (assuming the character of the particles in the jet 
remain constant) . Integrating the monitor signal over the total 
running time T of the jet we have: 
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Figure 3.- Typical "5-micron minusil" particles 


LOGNORMAL FIT TO PARTICLE FREQUENCY 



05.10.15.2.25.3 .4 .5 .6 .7 .75.8 .9 1.0 l.l 1.2 1.25 1.3 


PARTICLE DIAMETER D (MICRONS) 
Figure 4.- Minusil particle density function 
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The mass scattering coefficient for the sample is: 

k s (p,A,e M ) e e s (p,A,e M )/p(e M ,t) = w M (p,A,e M ,t)K(p,A,e M )/p (e M ,t) 

where p ( ©14) denotes the mean mass density of particles in the 
common volume at scattering angle 0 and K(P,A, 0 j^) is the system 
calibration constant. 
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The jet calibration procedure yields a relation between the 
particle mass density in the common volume and the overall mass 
flow rate M(t) in the jet: 


P (0 >t) = C ( 0 ) M (t ) . 


Then 


k s ( P,A,0 M ) 

V " K(P,A, 0 m ) X C(0 M 
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M(t) 


k (P,A ,0 ) 

K(P,A, 0 m ) X C(0 M )M 


where M is the total mass of particles expended during the run. 
Then the mass scattering coefficient for the monitor is: 

2 

k s (p,fl,e M ) = v x k(p,a,0 m )/mxc( 0 m ), ™. ster - 

Once the monitor scattering coefficient is known the divided 
signal can be reduced since: 

K sc (P,A,6,t) k s (P,A,8) K(P,A, 0 m ) - (0 _ t) 

W M (P,A,e M ,t) “ k s (P,A,6 M ) x K (P , A , 6 ) x p( 0 M ,t) 


or 
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K(P,A,0) C(6 M ) . 
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These equations are the basis for the data reduction. All of 
the constants needed are derived from calibration of the jet and 
the nephelometer system. 

The preceding analysis assumes that only single scattering 
occurs in the common volume during the measurements. Single 
scattering was proved experimentally by taking data at different 
mass flow rates and confirming a linear relation between signals 
observed and mass flow. Additionally, calculations based on the 
measured mass scattering coefficients and particle densities 
confirmed that the optical thickness through the jet was always 
small enough so that only single scattering was possible 
(t < 0.15). 


Derivation of Matrix Elements 


The Stokes vector of the light scattered from the sample 
volume at any angle 0 and accepted by the analyzer in the photo- 
meter can be written in matrix form as : 
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The symbols [Aij] , Oijl and [Sj_j] denote the 4x4 Mueller 
matrices describing the photometer analyzer, light source 
polarizer and the scattering medium. The elements of the 
scattering matrix are functions of the particle size, shape, and 
material, as well as the angle between the incident and 
scattered beams. Only the first Stokes parameter of the 
scattered light can be measured by our photometer. However, by 
making measurements with suitable combinations of polarizers and 
analyzers, each of the sixteen possible matrix elements can be 
obtained. In these measurements, only linear and right circular 
polarizers are used. The polarizers were H, V, D, d, and R. The 
usual sign convention is used here. If one looks along the light 
beam in the direction of propagation, the symbol H denotes 
polarization parallel to the scattering plane of the instrument 
and V denotes perpendicular polarization. The symbol D denotes 
polarization at 45 degrees to the horizontal with the electric 
vector pointing upward and to the right; d is orthogonal to D, 
and R denotes clockwise rotation about the direction of propa- 
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gation when looking in the direction of propagation. For the 
polarizer, the symbols indicate the polarization form transmitted, 
while for the analyzer, they denote the polarization form accepted. 
For these measurements, the following combinations were used: 

UU, HH, W, VH, HV, DD, dd, Dd, dD , RR, 

Rd, dR, DR, RD, HD, DH, HR & RH. 

In addition, consistency checks were made between various 
measured and computed values. For example, for an unpolarized 
source and receiver: 

I 

k s (U,U, 0 ) = j— (U , U , 0 ) = S i;L (m,x , 0 ) 
o 

while some matrix manipulations give 

s H ( m,x o , e) = j[k s (H,H,0)+k s (V,V,0)+k s (H,V,0)+k s (V,H,0)J 

Cross checks of this type were made whenever possible and were 
found to be satisfied within our experimental error. 

The measurements confirmed that the following relations 
held between the off-diagonal matrix elements for all scattering 
angles: S 12 = S 21 , S34 = -S43, and S13 = -S 31 = S14 = S41 = 

S 23 = -S32 = S24 = S42 = 0 . The scattering matrix then has six 
independent non-zero elements: the condition for a symmetric, 
isotropic scattering medium. 

V.- RESULTS AND DISCUSSION 

Our results are summarized in Figs. 6 through 9 . The 
matrix elements Sn(m,x o ,0) and S22 (* , x 0 , 0 ) are shown in Fig. 6 
as functions of scattering angle. For comparison, the corres- 
ponding matrix elements for spherical particles with the same 
size distribution are also shown. For spherical particles the 
two elements S^ and S22 are equal. The experimental points 
shown represent averages of at least six distinct measurements. 

The flags represent one standard deviation of the data about the 
mean values. The mass scattering coefficient for unpolarized 
light (matrix element Sii(m,x o , 0 ) shows excellent agreement with 
the theoretical predictions for scattering angles of less than 80 
degrees. For larger scattering angles, the theory disagrees with 
the measurements. At backscattering angles, there is a full 
order of magnitude difference. These results confirm our earlier 
measurements for unpolarized light (ref. 4 ) . The better agree- 
ment obtained here for scattering angles less than 80 degrees 
must be in part due to our improved knowledge of the particle 
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Figure 6.- Matrix elements S, , 


and S 


99 


vs 0 




sizes used in the experiment and our calibration procedure. 


The fact that the measured matrix elements Sn(m,x o , 0 ) do not 
equal the element S 22 (ra> x O/ 0 ) is not surprising. That particular 
relation holds only for spherical particles, and the experimental 
data merely reflect the fact that randomly oriented irregular 
discs do not resemble spheres in every respect. In fact, recalling 
the more fundamental amplitude scattering matrix from which the 
intensity scattering matrix is derived (ref. 6), we have: 
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const x 
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when the reference plane and subscripts are as described earlier. 
For spheres, the amplitude scattering functions A3 and A4 are 
equal to zero. Only then do the scattering matrix elements, 

Sll = S22 and S33 = S44. Our results merely confirm that the 
amplitude matrix elements A3 and A4 are not zero in this case. 

The better agreement with the theoretical curve exhibited by 
S22 is puzzling. We have no explanation for this fact at this 
time . 


The measured matrix elements Si2(m,x o ,0) and S2l(m,x o ,0) are 
shown in Fig. 7 as functions of scattering angle 0 . These 
elements were found to be equal within experimental accuracy. The 
corresponding matrix element for the spherical particle model is 
shown for comparison. The measured results show a preponderance 
of vertical scattering (negative values for matrix element S12) / 
while the spherical model predicts a preponderance of horizontal 
scattering. These measurements show the same general behavior 
as the measurements on atmospheric samples reported by Rozenberg 
and Gorchakov (ref. 12 ) . 

Fig. 8 shows the measured scattering matrix elements 
S33 (m,x 0 , 0 ) and S44(m,x o ,0) as functions of scattering angle 0 . 
Again, the fact that these elements are not equal implies, as in 
the case of the first two diagonal elements, that the amplitude 
functions A3, A4 are non-zero. The close agreement of the 
measured elements to the theoretical prediction is encouraging, 
since the element 833(10, x o , 0 ) shows, at least in the Mie theory, 
much more sensitivity to index of refraction than does Sn or S22* 

Finally, the off-diagonal matrix elements S34(m,x o ,0) = 

-S43 (m,x 0 , 0 ) is shown in Fig. 9 . The equality of these two matrix 
elements was confirmed within our experimental accuracy. Again, 
the theoretical values for spherical particles are shown for 
comparison. The agreement in this case is somewhat poorer than 
for the diagonal elements. 
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Figure 8.- Matrix elements and vs 9 
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It is tempting to try to generalize from these limited 
results. However, the data presented here represent measure- 
ments on one particle material and shape, fitting one size 
distribution and at one effective wavelength. Generalizations 
will require much more information. At present, measurements 
are being made at other wavelengths and the entire set of 
measurements will be repeated for the same particle material and 
shape but with larger sizes. Samples of the material are 
available with three larger size ranges. Measurements of the 
scattering matrix elements for other particle materials and 
shapes are also planned. 

There are several conclusions to be drawn from these 
results, however. First, that the unpolarized mass scattering 
coefficient for a polydisperse system of irregular but randomly 
oriented particles shows a remarkable similarity to the 
corresponding coefficient for spherical particles, but that 
predicting backscattering for such particles on the basis of the 
spherical model can lead to serious errors. Second, that the 
particle dispersion system used in these experiments generates 
a scattering medium whose matrix exhibits the properties of an 
isotropic symmetric medium. 
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